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Turbo-Charging DERIVE XM Version 2.57 

Several users have commented that the XM (extended memory) version of DERIVE ran significantly 
slower than the conventional memory version. 

Our test showed that XM version 2.56 ran about 30% to 40% slower on average (on some problems it 
ran twice as slow) as regular DERIVE. However, on big problems, XM was many times faster (e.g. on 
computing 4000!) since it did not have to thrash memory like regular DERIVE. 

I have spent several weeks hand-assembling some of XM's critical routines to give the ultimate in 
performance. Thus I am happy to report that the just released DERIVE XM version 2.57 runs as fast or 
faster than the regular DERIVE on small problems and much faster on big ones! 

Also new with XM version 2.57 is our ability to limit the amount of extended memory XM uses. This 
is done by means of an MS-DOS environment variable you can set before starting DERIVE XM. 

Users of multi-tasking environments as Windows or OS/2 may find this useful to prevent XM taking 
all memory resources. 

Aloha, 

A1 Rich, Soft Warehouse, Inc. 

Members of the DERIVE User Group are eligible for a free exchange of their DERIVE XM license to 
DERIVE XM version 2.57. Please send your original DERIVE XM diskette and a reference to this 
offer to one of the following addresses: 

For Non-European members: For European members: 

Soft Warehouse, Inc. Soft Warehouse Europe 

3660 Waialae Avenue, Suite 304 Softwarepark 

Honolulu, Hawaii 96816 A-4232 Hagenberg 

USA AUSTRIA 


DERIVE - BOOK - SHELF 


Es freut mich, Ihnen ein neues DERIVE Buch ankiindigen zu konnen, das von DUG-Mitgliedern ge- 
schrieben wurde: 

Mauve-Moos, Mathematik mit DERIVE, Ferd. Dummler, Bonn (mit Diskette Dummler Nr. 4588) 


2nd Krems Conference on Mathematics Education, September 27 - 30 


Ideas for Future Versions of DERIVE; Success and Failure in Mathematics; Number Theory with 
DERIVE; The Assessment of Mathematical Ability in the "Light” of DERIVE; Approximation 
Methods via DERIVE; From Harmony to Chaos; Deterministische und Stochastische Simulationen fur 
die Schule; Learning Visually; Using DERIVE to Explore Meaningful Applications in Calculus; Whi- 
te Box-Black Box-Principle; A Mathematical Model of a Firebreak using DERIVE; The World Sys- 
tem of Johannes Kepler in Stereo-Vision using DERIVE; An Experience in Algebra and Discrete 
Mathematics; Using DERIVE in the Calculus Classroom; Geometric Maximum and Minimum Prob- 
lems; Using CAS in Teaching Mathematical Modelling; The Understanding of the Concept of a De- 
rivative; The Use of Graphics Calculators and CAS; The Application of a CAS as a Tool in College 
Algebra; Symbohc-Computation-unterstutzter Mathematikunterricht; The DERIVE Project in 
Petropolis; DERIVE-centred Research at the University of Plymouth; Methodische und didaktische 
Bemerkungen zum Einsatz von DERIVE; DERIVE im Unterricht der 11. und 12. Schulstufe. 








D-N-L#ll 


LET TEH 


O f 


THE TV iron 


p 1 


Liebe DERIVE Anwender! 

Viele von Ihnen haben zum Teil recht aus- 
fuhrliche Ruckmeldungen auf meine Fragen zur 
Gestaltung des DNL und der Arbeit der DUG 
zuriick geschickt. Dafiir mochte ich mich 
nochmals recht herzlich bedanken. Ich habe nun 
die Anregungen, die an mich herangetragen 
worden sind, in einer Wunschliste zusammen- 
gestellt. Einige Themen wurden mehrfach ge- 
nannt. Die Liste ist ungeordnet. Manches wurde 
bereits in diesem Jahr und einiges wird in die- 
sem neuen DNL behandelt. Einiges ist in Vor- 
bereitung (Abbildungsgeometrie mit Matrizen, 
Arbeitsblatter, Anwendungen fur die Elektro- 
technik - H. Scheuermann, ...), aber viele inte- 
ressante Ideen bleiben offen! Ich danke auf 
diesem Weg fur die vielen zustimmenden Zu- 
schriften und fur ihre Beitrage. Auch in diesem 
Heft findet sich eine recht internationale Mi- 
schung: USA, UK, Portugal, Slowakei und die 
BRD. 

Ende September findet in Krems, Niederoster- 
reich, die 2. DERIVE Konferenz statt. Auf der 
linken Seite konnen Sie die fur diese Konferenz 
angenommenen Vortrage finden. Wir sind si- 
cher, dass auch von diesem Treffen wieder 
wertvolle Anregungen und Impulse ausgehen 
werden. In den nachsten DNLs werde ich dar- 
iiber berichten. 


Dear DERIVE Users, 

Many of you have answered my questions 
concerning the development of the DNL and 
the further work of the DUG. Many thanks for 
that. I collected all your ideas and suggestions 
in a wishlist. There is no special order in this 
list. Some items were mentioned several times. 
Some contributions of 1993 are meeting your 
wishes, some answers are in preparation (Map- 
pings with matrices, worksheets, applications 
for electrical engineering - H. Scheuermann, 
...). But many interesting ideas are still open. I 
am very indebted for so many positive letters 
and for your contributions. I think that we have 
a very international mixture in this issue, too: 
USA, UK, Portugal, Slovakia and Germany. 

By end of September the 2 nd DERIVE Confer- 
ence will be held in Krems, Lower Austria. We 
are sure that this meeting will bring many inspi- 
rations and impulses. I will report in the next 
DNLs. 

On the left page you can find the lectures which 
have been accepted for this conference. 

with my best regards 


Mit den besten GriiBen Ihr 


DERIVE USER GROUP Members’ Wishlist 

Teaching Algebra, Calculus, ODEs, Integration; Vektoranalysis; Economics; Sta- 
tistics; Probability Theory; Okonometrie; Good examples which require CAS for 
their solution; Astronomical Content; Celestial Mechanics; Chaos; Technical, 
Physical and Economic Applications; Approximation; Ingenieurwissenschaften; 
Advanced Stuff; Worksheets; Matrizen; Abbildungsgeometrie mit Matrizen; 2D- 
und 3D-Plots; Anwendungen zur analytischen Geometrie; Anregungen zu Printer 
& Plotter; Unterrichtskonzepte - Beispiele; Erfahrungsberichte aus der Schule; 
Erfahrungen liber den Gebrauch von DERIVE bei wissenschaftlichen Arbeiten; 
Technische Problemlosungen (zB Elektrotechnik); Didactics; 3D-Graphing of 
Complex Numbers; Mathematical Games; Exporting Plots into other Programs; 
Recursive Programming; Differentialgleichungen - Differenzenverfahren; Back- 
ground Knowledge to Implementations of DERIVE 
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The Derive-News-Letter is the 
Bulletin of the Derive-User 
Group. It is published at least 
three times a year with a con- 
tents of 30 pages minimum. The 
goals of the D-N-L are to en- 
able the exchange of experi- 
ences made with Derive as well 
as to create a group of experts 
to discuss the possibilities of 
new methodical and didactic 
manners in teaching Mathemat- 
ics . 


Editor : Mag. Josef Bohm 
A-3042 Wurmla 
D ' Lust 1 
Austria 

Phone/FAX: 43- ( 0 ) 2275/8207 
e-mail : no j o . boehmOpgv . at 


Contributions : 

Please send all contributions 
to the below address. Non- 
English speakers are encouraged 
to write their contributions in 
English to reinforce the inter- 
national touch of the D-N-L. It 
must be said, though, that non- 
English articles will be warmly 
welcomed nonetheless. Your con- 
tributions will be edited but 
not assessed. By submitting ar- 
ticles the author gives his 
consent for reprinting it in 
the D-N-L. The more contribu- 
tions you will send to the Edi- 
tor, the more lively and richer 
in contents the Derive-News- 
Letter will be. 


Preview: (Contributions for the next issues) 

Fluid Flow in DERIVE, Reuther a.o., BRA 

Newton's Chaos, Graphic Integration, J. Bohm, AUT 

Computer Aided Mathematics in School, Karl-Heinz Keunecke, GER 

DERIVE in Hawaiian Classrooms, E. Sawada, USA 

Applications in Electrical Engineering, H. Scheuermann, GER 

Stability of Systems of ODEs, A. Kozubik, Slovakia 

Minimization of a „Flat Function 44 , C. Lopes, Portugal 

Maximum and Minimum Values - A Tool for Teachers, E. Zott, AUT 

(DNL#12 will be published December 1993) 


Impressum : 

Medieninhaber : DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrif t 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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The function MAP may be changed in an obvious way if the vector v is an array of higher 
dimension. 

You cannot use this MAP-function from 1993 now because there is a function MAP imple- 
mented in the recent DERIVE versions. Functions map_list (u, x, c) and 
map list (u, k,m, n, s) are very useful - I recommend a glimpse into the Online-Help. 
map (u, x, c) returns the truth value true. (Josef 2006) 


#1: MAP_LXST(LN(x) s x, VECTOR (k, k, 10)) 

*2: [0, LN(2), LN(3), 2*LN(2), LN(5), LN(6) , LN(7), 3*LN(2), 2*LN(3), LN(IO)] 

#3: MAP(LN(x), x, VECTOR (k, k, 10)) 

# 4 : 

Glynn D Williams, Gwynedd, UK 

I regularly use the DERIVE program, but find it annoying that a different set of keys is needed to edit 
an algebraic expression than to select the expression or subexpression to be edited. The reason given 
in the manual for these inconsistencies is that certain combinations of keys which would be logical are 
not recognized; e.g. — , HOME and END for editing an expression in the bottom editor line; or 
CTRL+t and CTRL+>1 for moving the cross in the plot window one large division. However, if 
ANSI.SYS. or a compatible public-domain or shareware driver is loaded, it is quite easy to patch the 
keyboard so that it behaves logically. The batch file below does this: 


esc 

0; 71; 17; 19p 

rem 

Home to 

give A Q A S 

esc 

0 ; 79; 17 ; 4p 

rem 

End to give A Q A D 

esc 

0; 75; 19p 

rem 

Left to 

give A S 

esc 

0 ; 77 ; 4p 

rem 

Right to 

give A D 

esc 

0; 141; 0; 73p 

rem 

Ctrl-Up 

to give PgUp 

esc 

0; 145; 0; 81p 

rem 

Ctr-Down 

to give PgDn 

der. 

ive . exe %1 




esc 

0 ; 7 1 ; 0 ; 7 lp 




esc 

0 ; 7 9 ; 0 ; 7 9p 




esc 

0 ; 7 5 ; 0 ; 7 5p 




esc 

0; 77; 0; 77p 




esc 

0; 141; 0; 141p 




esc 

0; 145; 0; 145p 





Six keys need to be patched: they are all in the extended code set, and they are all to be 
patched to give codes in the normal code set. Note that they must all be "unpatched" after exiting from 
DERIVE, as otherwise bizarre results may be caused especially if an enhanced command interpreter 
such as 4DOS or Doskey (or a word-processor or editor) which use the arrow keys is in use. 

esc refers to a small machine-code program esc . com, which outputs the ESC character followed by 
a [ to standard output when it is activated., and then tags on the characters given as parameters. With 
this patch in operation one loses the ability to select a sub-expression of a highlighted expression 
while there is an expression in the edit window, but this is rarely needed; apart from this there is no 
ambiguity. If it is needed one can always isolate the sub-expression in a spare numbered ’’slot” (by 
highlighting it and then copying it using F3) before editing it; there is usually enough memory to do 
this. 
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DNL: In DERIVE 2.55 the F 6-key enables comfortable editing now , but I am sure that there are a lot 
of our members who don 7 own the latest version. I find it very interesting to see how to adapt a pro- 
gram for one's personal use. 


A.C. Robin, Colchester, UK 

I had some difficulties trying to use DERIVE for solving some problems, and I wrote to Soft 
Warehouse about them. I am enclosing copies of the letters I sent together with a reply I had to the 
first difficulty. I thought that you or readers of the DNL might be interested: 

'0 1 ' 

1 2 


+ c *x + d*x + ej , ma4) 

2 

2220 * x 7 699 * x 35243 

+ + + 

4S49 13187 35244 

See David Stoutemyer's answer: 

The problem disappears if before simplifying the FIT example you choose 

Options Precision 15 

The explanation is that FIT uses the straightforward ’’normal equation” method together with 
appproximate mode, and this mode typically loses about half the working precision, which defaults to 
6 digits. Also, using the monomial x m as basis functions is notoriously subject to additional roundoff 
error. 

Methods such as singular- value decomposition together with orthogonal polynomial bases are 
more accurate, but they require an impractical amount of code for a general purpose computer algebra 
system intended to run in as little as 512 kilobytes. Thus, given the arbitrary precision of DERIVE, it 
is more appropriate merely to increase the working precision as the desired polynomial degree in- 
creases. 

Sincerely yours, ... 

4 3 2 

x x ll*x 7*x 

# 4 : + + + 1 

24 12 24 12 

This is no problem with recent versions of DERIVE: 

#5: PrecisionD-ig-its = 3 

#6: NotationDigits = 3 

4 3 2 

[432 1 xx 11.x 7.x 

#7: FrT(Lx, a*x + b*x + c*x + d*x + ej , ma4) = - + + + 1 

24 12 24 12 


The FIT-function does not give back the 
correct quartic. 


# 1 : 


#3: 


ma4 = 


#2: FIT 


2 4 

3 8 

4 16 

L 5 31 J 

[ * 

(Lx, a*x + b. 


338 . x 


8115 


787.x 


9462 
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Mr Robin’s 2nd problem: 

I wonder why DERIVE 2.55 gives the following, have I done something wrong or is there a bug? 

As you can see in the output below (left column, editor), I defined the cube function as x A 3, 
and defined the d cube -function as its derivative, namely 3x A 2. When I asked DERIVE to integrate 
x* dcube (x) and (x+u) * dcube (x) , it gave sensible and correct answers (u was assumed con- 
stant as I wished). However when I asked it to find the integral of x* dcube (x+u) then a strange 
expression occurred, and this was also obtained when just dcube (x+u) was simplified. DERIVE 
does correctly simplify dcube ( z ) as 3z A 2, so why does not dcube (x+u) simplify to 
3 (x+u) A 2? 

I came across this problem by trying to use a much more complicated function than x A 3, so 
that I really needed to define the function and its derivatives in this way. 


3 

#1: CUBE(x) := x 

d 

#2: DCUBE(x) := — CUBE(x) 

dx 

#3: J x • DCUBE(x) dx 

4 

3 • x 

#4: 

4 

#5: J (x + u) • DCUBE(x) dx 

3 

x • (3-x + 4 • u) 


In 1993 I overcame the problem introducing an aux- 
iliary variable. In the meanwhile this is not neces- 
sary, because the later versions of DERIve behave 
as expected, Josef 

3 

#1: CUBE(x) x 

d 

#2: DCUBE(x) := — CUBE(x) 

dx 

4 

3 -x 

#3: J x • DCUBE(x) dx = 

4 


4 

3 • x 3 


#7: 

J x«DCUBE(x + u) dx 



d 

#8: 


x- CUBE(x + u) dx 


• 

d x + u 

#9: 

DCUBE(z) 



2 

#10: 

B-z 

#11: 

DCUBE(x + y) 


d 


#4: J (x + u) • DCUBE(x) dx = 

4 

2 

#5: DCUBE(x + u) = 3-(x + u) 

2 2 
#6: 3-(x + u) -x = 3-x-(x + u) 

3 

f 2 (x+u) 

#7: J 3-x-(x + u) dx = 


+ u *x 


.(3.x - u) 
4 
3 


#12: CUBE(x + y) 


(x + u) .(3*x - u) 


d x + y 


# 8 : 


J x • DCUBE(x + u) dx = 


4 
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Message 3000 was entered on 5/21/93 at 2:06 PM 

From SOFT WAREHOUSE to PUBLIC about CHAINE RULE DIFFERENTIATION 


Recently there has been much discussion on this BBS about having DERIVE use 
the chain rule when simplifying partial derivatives of arbitrary functions. 
The most general case is derivatives of the form 

- 7 — F (G (x,y) ,H (x,y) ) . 
dx 

Applying the chain rule to this yields an expression of the form 


r 


v 


lim 

01— >G (x,y) q[0 ]_ 


\ 

F ( 01 , H (x,y) ) 

) 


d 

G (x,y) + 

dx 


f 

V 


lim 

@2— >H (x,y) (j_0 2 


\ 

F (G (x,y) , 02 ) 

J 


d 

H (x,y) . 

dx 


Of course if G and H are specified rather than arbitrary functions, simpler 
results would be produced. For example, DIF (F ( sin (x) ) , x) would simplify to 
LIM (DIF (F (01) , 01) , 01, SIN (X) ) *COS (X) . 

If DERIVE performed such transformations, would this help satisfy your dif- 
ferentiation needs? 


Message 3021 was entered on 5/28/93 at 2:06 AM 

From ROGER FOLSOM to SOFT WAREHOUSE about DERIVE SUGGESTIONS 

CHAIN RULE DIFFERENTIATION. You ask whether it would meet user needs if De- 
rive calculated 

F (G (x,y) ,H (x,y ) ) as ... (Please see message 3000.) 
dx 

I think this would be useful for those working with specified functions. 

But since I work mostly with what Derive calls arbitrary, unspecified, 
functions (which I'm used to calling general functions - which phrasing, 
arbitrary or general, is more commonly used among mathematicians?) , I fear 
that carrying around all those limit operations would chew up an enourmous 
amount of memory. Consequently, I think I'm better off defining F (01, 02) :=, 
G(x,y) :=, H(x,y) :=, and for the system 

G (x, y) - 01 = 0 
H (x, y ) - 02 = 0 

setting J : = JACOB IAN ( [G(x,y)-01,H(x,y)-02] , [01,02,x,y] ) , noting that 
J . [d01 , d02 , dx, dy] = [ 0 , 0 ] , and if JN is the first two columns of J and JX 

is the last two columns of J, then solving 
JN . [d01,d02] = - JX [dx,dy] 

gives [d01,d02] = - JN A (-1) . JX [dx,dy] 

Substituting this last result into the total differential of F (01, 02) that 
is into GRAD (F ( 0 1 , 02 ) [01,02]) . [d01,d02], gives the partial derivative of 

F(.) wrt x (if dy=0) or y (if dx=0) . 

But it would be wonderful if DERIVE automated this process. If 01 = G(x,y) 
and 02 = H(x,y), then the objective should be to end up with something at 
least as simple notationally as: 

d d dH d dH 

F(G(x,y) ,H(x,y) ) = F (01,02)* + F (01,01)* 

dx d 01 dx d 02 dx 

Note that much memory would be saved (I think) and certainly much screen 
and paper real estate would be saved if DERIVE were willing to substitute 
the symbols F for F(01,02) and H for H(x,y) . 
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Message 3034 was entered on 6/2/93 at 10:29 PM 

From SOFT WAREHOUSE to ROGER FOLSOM about #3021/DERIVE SUGGESTIONS 

CHAIN RULE DIFFERENTIATION. Your notation for the partial derivative of 
F(G(x,y) ,H(x,y) ) wrt x appears compact only because of the hidden assump- 
tions it makes. Specif ically, the assignments 

@1 := G(x,y) and @2 := H(x,y) 

would have to be made as well as the declaration that G and H are functions 
of x. Note that there are a couple of typos in your example: the first 
dH/dx should be dG/dx and the last @1 should be @2. 

Aloha, A1 Rich, Soft Warehouse, Inc. 


Message 3045 was entered on 6/6/93 at 4:41 PM 

From ROGER FOLSOM to SOFT WAREHOUSE about DERIVE SUGGESTIONS 

CHAIN RULE DIFFERENTIATION: you are quite right that in my notation for the 
partial derivative of F (G (x, y) , H (x, y) ) wrt x, namely 

d d dG d dH 

F(G(x,y) ,H(x,y) ) = F(X,Y)* + F(X,Y)* 

dx dX dx dY dx 

(where I hope I have fixed the typos in my earlier message), with F(X,Y) :=, 
G(X,Y) : =, H(x,y) : =, I "forgot" to state the assignments 

X := G(x,y) and Y := H(x,y) 

Prior to today, my previous experience was that if I had included those 
last two assignments. Derive (2.08) would have refused to differentiate 
F(.) with respect to X and Y, which is why I simply kept those assignments 
"in my head". Today, however on checking this issue, I found that Derive 
was willing to differentiate with respect to X, Y, x and y, so I must have 
had some confusion in my earlier work (possibly an uncleared conflicting 
assignment) . Or maybe I'm confused now. 

Aside fom that issue, note that even if the approach I suggested requires 
as many characters or keystrokes as the limit notation you suggested, my 
result line is shorter (even if I need more lines for function assign- 
ments) , and therefore easier to read, print, and manipulate. 

(Your discussion used @1 for X and @2 for Y. Does @ have some special mean- 
ing? I don't find it in the manual index, or in the manual itself.) 

Message 3000 was entered on 5/21/93 at 2:11 PM 

From GREG SMITH to JACKNOL about #2998 / MATH TOOLKITS FROM U OF AZ 

We offer some of the Arizona toolkits for download here as a service to us- 
ers of the BBS, but it would be too labor-intensive to offer them on disk. 
You may want to download one to get the address and then write or call the 
authors to see if they are available some other way. 

Message 3014 was entered on 5/2693 at 11:17 AM 

From DELAWARE to PUBLIC about UNIV. OF ARIZONA PROGRAMS 

To all those interested in the U of AZ software: I recently bought them 
all, meaning I sent the UA Math Dept $28.00 for 14 diskettes which include 
both the toolkits (9 diskettes) and non-toolkits (5 diskettes) . The latter 
include slide shows, simulations. Are You Ready quizzes, etc. I find many 
of them very useful, both as assignments and in class preparations. De- 
scriptions come with the diskettes, in fair detail, but the on-line in- 
structions are usually insufficient. And, since this is freeware, I've 
passed them around my region (KC MO) . I recommend them all. Richard Dela- 


ware 
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Message 3036 was entered on 6/3/93 at 10:16 PM 

From DELAWARE to PUBLIC about UNIV. OF AZ FREEWARE 

To those interested in obtaining the UA freeware, the info is: 

Dept of Mathematics, Univ. of Arizone, Tucson, AZ 85721 

Phone: (602) 621-6893. Call and have them send you an information sheet 

which has all the cost, options, etc. data. Goodluc. RD 

Message 3013 was entered on 5/26/93 at 7:09 AM 

From JERRY GLYNN to BOOM-BOOM about S A N/ (S A N+C A N) 

I sent a note to A1 about integral ( sin (x) A n/ ( sin (x) A n+cos (x) A n) , x, 0 , tc/ 2 ) . I 
noticed that it seemed to be = tt/ 4 for all n>0 so how about if Derive would 
recognize this ... he said that would be too narrow for a separate rule. By 
the way Derive does better than Mathematica on the problem. Why is the an- 
swer 7i/4? It sounds like it's the area of a circle who's radius is 1/2 but 
how does that connect? 

Message 3023 was entered on 5/30/93 at 12:10 PM 

From HADUD to JERRY GLYNN about COMMENT ON #3013 

The integral in question illustrates what happens when you ignore the spe- 
cial symmetry properties of a problem and try to solve it with a general 
program like DERIVE. In the present case the integrand is complementary, 

i . e . 

f (x) + f (7i/2-x) = 1 

Integrating this equation and substituting u for (7i/2-x) in the second in- 
tegral you find immediately 

2 INT (f (x) ,x, 0 ,tu/2) = n/2. 

This shows, incidentally, that the result holds for any finite n (positive, 
negative real, complex) . DERIVE could in principle, test for complementar- 
ity and thus solve problems of this kind easily. Most people would agree, 
though, that this is too special a case to build into a general program. 

The lesson to be learned here is that even a great symbolic program like 
DERIVE is no substitute for mathematical insight! 

Message 3024 was entered on 5/30/93 at 1:12 PM 

From JERRY GLYNN to HADUD about #3023 / COMMENT ON #3013 

Thanks very much for your comment. Where do I loo for more examples of 
mathematical insight? Can we do a few more examples here? 

Message 3031 was entered on 6/2/93 at 10:02 PM 

From HADUD to JERRY GLYNN about MORE ON COMPL . INTEGRANDS 

Pursuant to our telephone conversation here are some more thoughts on 
evaluating definite integrals with complementary integrands. The general 
identity applicable in these cases is 

INT ( f (x) , x, a, b) = INT ( 1/2 ( ( f (x) +f (a+b-x) , x, a, b) (1) 

This holds for any f (x) but is useful only if the modified integrand on the 
right is simpler (i.e. integrable) than the original f (x) . This is the case 
in your example with f (x) =sin (x) A n/ ( sin (x) A n+cos (x) A n) , a=0 b=7u/2. To apply 
this under DERIVE use the right-hand side of (1) to generate an expression 
(after declaring f (x) :=) . Then set f (x) equal to your function and simplify 
the expression with proper values for 'a' and 'b' . Further examples (all of 
which DERIVE gives up on) are: 










1 
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Linear : 
Reciprocal : 


g(x) = a + b - x p = 1/2 

g (x) = (a b) /x p = 1/2 

special case: a -> 0, b -> inf. 


(a b) 


INT ( f (x) , x, a, b) = INT ( 1/2 ( f (x) +f (a+b-x) ) , x, a, b) 

INT ( f (x) , x, a, b) = INT ( 1/2 ( f (x) + (a b) /x A 2 f (a b) /x) ) , x, a, b) 


3) IMPLEMENTATION in DERIVE: The functions INTI (f (x) , x, a, b) and 

INT2 ( f (x) , x, a, b) implement the above cases. Their syntax is the same as 
for INT. Thus if INT fails, before slashing your wrist try INTI and then 
INT2 . 

Function definitions: 

fF(x):= 

INTl(fct,x,a,b) ELEMENT([fF(x):=fct , J(l/2(fF(x)+fF(a+b-x)) ,x,a,b)] ,2) 
INT2(fct,x,a,b) 

ELEMENT([fF(x):=fct,IF(a=0,IF(b=oo,/(l/2(fF(x)+l/x A 2fF(l/x)),0,oo) J(fF(x), x, 0, b) , 

J (l/2(fF(x)+l/x A 2f F(l/x)) , x, 0 ,»)) , J (1/2 • (f F(x)+a- b/x A 2f F(a- b/x)) , x , a, b) , 
J(l/2(fF(x)+a*b/x A 2fF(a*b/x)) ,x,a,b))] ,2) 

Comments: a) Since you cannot have a function as an argument in a 

DERIVE function definition, the above definitions use 2- 
element vectors. The 1st element assigns the given inte- 
grand expression (fct) to an auxiliary function fF (x) . 

The 2nd element then performs the integration and is re- 
turned as the result. 

b) The definition of INT2 contains the special case of a=0 
and b=inf, hence the IF functions. 

(See the comments of 2006 at the end of this message, Josef) 


4) EXAMPLES (in all these cases INT fails) : 

1 . INTI (SIN (x) A n/ (SIN (x) A n+cos (x) A n) , x, 0 , pi/2 ) 

2. INTI (LN(u) /LN(u (u-1) ) ,U, 0, 1) 

3. INTI ( (3x A 2-3x+l) LN (x) /LN (x (1-x) ) , x, 0, 1 ) 

4. INTI (x SIN (x) / (2+SIN (x) ) ,x, 0 , pi) 

5. INTI (SQRT (l-x A 2) / (1+EXP (q x)),x,-l,l) 

6. INTI ( ( 1+x) A 3 COS (x) / (l+3x A 2) ,x, -1, 1) 

7. INTI (EXP (-x A 2) SIN(x) ,x,-c,c) 

8. INT2 (1/ ( (l+x A 2) (l+x A n) ) ,x, 0, inf) 

9. INT2 (LN (x) / (l+x A 2) ,x, 0, inf) 

10. INT2 (x ATAN (x) / ( (l+x A 2) (l+x+x A 2) ) , x, 0, inf ) 


Comment: All the above cases are contrived examples and of little prac- 

tical use. However INTI does solve the general case where f (x) 
is odd and a=-b . DERIVE 2.55 does not (in general) recognize 
that this gives a zero result (see example 7.) . This is a se- 
rious shortcoming and makes a strong case for incorporating 
INTI into the INT function. 
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Comments of 2006: 

First of all: Recent versions of DERIVE don't have problems with most of the integrals men- 
tioned above. DERIVE has learnt a lot between version 2.55 and later ones. Moreover using 
the wonderful "Stepwise Simplification" tool you can find the rules implemented as suggested 
in this email: 


tt/2 


#14: 


SIN(x) 


n n 

SIN(x) + COS(x) 


dx 


/ (F(x) + F(a + b - x)) dx 

b a 

/ Ffx) dx 

a 2 


n/2 


#15: 


#16: 


1 

— dx 
2 


8 


#17: 


- x 

c *SIN(x) dx 


-c 


/ (F(x) + F(a + b - x)) dx 

b a 

/ F(x) dx 

a 2 


#18: j 0 dx 


-c 


#19: 0 


Trying the examples using the multiplicative equivalent, we can still observe problems: 

<L> 

LN(x) 


# 20 : 


1 + x 


dx = 0 


# 21 : 


2 n 

(1 + x ).(1 + x ) 


dx = 


n 2 2 

X * (x + 1) + X + 1 


dx 


2*«/3 * 


ATAN 


V3*(2*x + 1) 1 


dx 


#22: 

<x> 

x • ATAN(x) 

dx - 

0 

x + 1 

2 

' 1 



2 2 
(l + x)*(i + x + x) 


3 


. s 

6 > 


Numeric integration of #22 gives: 


*24: 0.2839974233 
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HADUD tried something like programming with DERIVE 2.55, which worked in 1993. Recent 
version don't accept this way of assigning auxiliary variables in the first element of a vector 
which should be processed in the second one. We could write a small program or we avoid 
the auxiliary variable, substituting for x in fct (without using SUBST, to remain in the style of 
version 2.55). 


First attempt with the original INTI -function 


n/2 


*3: 


INTI 


n 



fF 

' n 

v 2 

n/2 

dx + J fF(x) dx 

0 

SIN(x) 

JT 

v n 


0 




l SIN(x) + COS(x) 

I rewrite this function and use it twice: 


#4: INTI (fct, x, a, b) = 


. (fct + lim fct) dx 

2 x-a + b - x 


*5 : INTI 


n 

SIN(x) n 

, x, 0, 

n n 2 


l SIN(x) + COS(x) 


n 

4 


*6: INTI 


(1 + x) *COS(x) 


1 + 3-x 


x, -1, I 


= 2*SIN(l) 


*7: 


(1 + x) .COS(x) 


1 + 3*x 


dx = 2*SIN(l) 


-1 


Same happens with INT2 

INT2 from 1993 fails: 


# 9 : 


INT2 


LN(x) 


x, 0, « 


2 1 ^ 

x *fF(x) + fF 


V A J 


dx 


1 + x 
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I rewrite this function, too and try once more and then I try the two Integrals which DERIVE refused to evaluate: 


IMT2(fct, x, a, b) = 
If a = 0 


If b = « 


#10: 

/ (1/2 • (fct + l/x A 2.LH1(fct, x ; 
J(fct, x, 0, b) 

. 1/x)), 

x, 0, 

«) 


J (1/2* (fct + l/x A 2.LB1(fct, x ; 

i 1/x)), 

x, 0, 

«) 


/(l/2*(fct + a*b/x A 2*LU1(fct, x, 

a*b/x)). 

. x, a. 

. b) 


J(l/2*(fct + a*b/x A 2*LU1(fct, x, 

a*b/x)). 

i x, a. 

, b) 


INT2 


* 11 : 


LN(x) 


1 + x 


x, 0, « 


= 0 


INT2 


* 12 : 


x 2 ).(l 


x, 0, 


n 

X ) 


n 

4 


#13: 


1 n 

dx = — 

10 4 


(1 + x )♦(! + x ) 


#14: 


2 11 

(1 + x ).(1 + x ) 


dx = 


2 11 

(x + lMx + 1) 


dx 


#15: 


#16: 


INT2 


INT2 


x*ATAN(x) 


2 2 
^(l+x)*(l+x+x) 


X, 0, *> 


2 

= JT * 


73 


IS IS J 


x*ATAN(x) 


2 2 
^(l+x)*(l+x+x) 


x, 0, *> 


= 0 . 2S3997423S 


It is interesting that DERIVE integrates #13, but does not accept exponents >10 in the sec- 
ond factor of the denominator. INT2 resolves the general case. Approximating the result of 
#15 and comparing it with the result in the last line on page 12, we can see that HADUD's 
functions are still working considering the adapted syntax. 
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(These lines can be saved in a DERIVE file for future use. At that time, they can be loaded into 
DERIVE via the Load command.) Collectively lines #1 through #6 define a 4 x 3 matrix which identi- 
fies the interval [ a , b\, its midpoint, m, and the values of the function, p(x), at a , m and b. Once p(x) 
and the values of a and b have been defined, line #6, when Simplified, will display the data matrix. 

To begin the bisection process, define the function, p(x), and choose values for a and b. For ex- 
ample, one can approximate V 3 by applying the algorithmus to the function p(x) = x - 3. Since 
p( 1) = -2 and p( 2) = 1, start with a = 1 and b = 2. Append the following line to those above. (Again 
DERIVE will provide the line numbers.) 

#7: [P(x) := x a 2 - 3, a := 1 , b := 2] 

Now Simplify line #6 (use the arrow keys to move the lightbar to line #6 and then execute the 
Simplify command). The display will be line #8 (see figure 1.) 

Notice that p(x) is negative at x = 3/2, the midpoint of [1 , 2], and also at a = 1 . Thus, the value 
of the left endpoint, a , should be reassigned to 3/2. That is, by changing the value of a to 3/2, the next 
interval \a 2 , b 2 ] = [3/2, 2] is defined. The procedure is then repeated, a new matrix is generated and the 
interval is subsequently refined. Author a new line (#9): 

#9: a := 3/2 

and then Simplify line #6 again. The new matrix, line #10, is shown in figure 2. This time note that the 
new midpoint, m 2 , is 7/4, and that the sign of /?(7/4) agrees with that of the right endpoint, b = 2. To 
reassign the value of b to 7/4, Author line #11: 

#11: b := 7/4 

and then Simplify line #6 again. Shown in figure 3 is the matrix corresponding to the interval, 

[3/2, 7/4]. 

Line #6, will by this time, have risen out of the viewing screen. It can of course be recovered by 
using the arrow keys to move the lightbar up to line #6, but it is usually more convenient to have it in 
the visible screen. To do this move the lightbar to line #6, press "a" for Author and then the F3 key 
followed by "RETURN”. (Pressing F3 does now the full job.) This procedure "pulls down" the 
highlighted line into the visible screen with a new line number. Now a matrix of values will be dis- 
played whenever the new line is Simplified. 

To display the data in decimal, rather than fractional form, move the lightbar to the most re- 
cently displayed matrix and press "x" for Approximate, follwed by "RETURN". The approximated 
matrix , line #13, is figure 4. (Press the ^-button or Ctrl+G in recent DERIVE versions. Ctrl+G 
has the additional advantage of setting the accuracy.) 

If you prefer to see only decimals as sequential matrices are generated, use the Options Preci- 
sion command to set the mode to "approximate" and set digits to the desired level. Shown here are six 
digits, the DERIVE default. While in approximate mode, the next matrix, that for [a, b ] = 

= [1.625, 1.75], is shown in figure 5. 

Although, in the example, the process will never arrive at the exact value of V 3, the midpoint m p 
will eventually provide a very good approximation. Note that the length of the interval containing the 
root is halved at each stage. Therefore, if the midpoint, mj, is used to approximate the root of p(x) after 
the / h step, the error of approximation will not exceed ( b\ - a\)!2 j . In the example, a x = 1 and 
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b\ = 2, thus the midpoint mj, is never further than 1/2 7 from V3. From the next matrix, we have 
[a 4 , Z? 4 ] = [1.625, 1.75] and thus m 4 = 1.6875, approximates V 3 with a maximum error of 2' 4 = 0,0625. 

To approximate a root to some predetermined accuracy, 8, find the first integer, n , such that (b\- 
a\)l2 n < 8 and then determine the midpoint m n , of \a m b n ] which provides the desired approximation. In 
the example, if we set 8 = 0.001 and note that 2~ n < 0.001 for n = 10, then m\o = 1.73241 approximates 
V3 with the desired accuracy, 0.001. 



a 

m 

b ' 


a 

m 

b ' 



a 

m 


b ' 



3 



3 

7 




3 

13 


7 


1 

- 

2 


— 

— 

2 



- 

— 


- 



2 



2 

4 




2 

8 


4 

p(a) 

p(m) 

p(b) 


p(a) 

p(m) 

p(b) 



p(a) 

p(m) 

P(b) 


9 

3 



3 

1 




3 

23 


1 

. 

~ L. 

4 

_L 


4 

16 

_L 



4 

64 

16 J 



fig 1 




fig 2 





fig 3 




- 

a 

m 

b ' 



■ 

a 


m 

b 

- 




1.5 

1.625 

1.75 



1.1 

625 


L . 6875 

1.75 




p(a) 

p(m) 

p(b) 



p(a) 


P Cm) 

p(b) 




L- 

0.75 

-0.359375 0.0625 J 



. -0.359375 

-0.152343 

0.0625 




fig 4 fig 5 

The DERIVE Algorithm in the Classroom 


The DERIVE bisection algorithm described here requires little practice to be mastered. It is ide- 
ally suited to classroom demonstration. After a short chalkboard introduction to the method, it be- 
comes clear that "hand by hand” calculations are prohibitive. At that point, engage the class in an in- 
teractive session with the microcomputer. At each stage of the process, have the class decide how to 
refine the interval which will be used to generate the next matrix. Later students can apply the method 
to other problems in the microcomputer lab. They will gain insight into the method, both how and why 
it works, as well as gain experience with DERIVE. There are a number of other opportunities which 
are made available by the algorithm. For example, have the students who understand the procedure 
and who are sufficiently experienced with DERIVE dissect and interpret the lines to compute the ma- 
trix at each stage of the process. Have them explore possible alternatives to the coding and possibly 
discover a more informative matrix. Ask more advanced students to research the version of the Heine - 
Borel Theorem which guarantees that an infinite nested sequence of closed intervals will have a non- 
empty intersection, {x*}. Have them argue, if not prove, that the value of x* must actually satisfy, pipe*) 
= 0, and that mj —> x* as j —> oo. Challenge them to provide examples which demonstrate the necessity 
for the nested intervals to be both closed and bounded. 

In general there are many opportunities which are provided by the method, especially in light of 
the simplicity with which DERIVE casts the calculations and ease of display. 


Fulks, W (1969). Advanced Calculus: An Introduction to Analysis (2nd ed.). New York: Wiley 
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Another challenge could be to enclose the whole procedure in a “program”. Without introduc- 
ing the programming features of Derive one can emphasize the recursive nature of the 
method by applying the ITERATES-command which cannot be offered too often to the stu- 
dents. I present one first approach including the given accuracy in the last column. This tool 
does not take into account any special case mentioned above (eg. f(x)=0 at one of the three 
positions) 


bisect (fct, 



f 

f 

f 


a, b, n) := ADJOIN 

[a, m. b, max_error], ITERATES 

IF 

SIGN( 1 im fct) 

= SIGN( lim fct), 


V 

V 

^ x-vjl 

x*v.i2 


v + v v 

2 3 4 

v , 1 v , 

.2 2 3 2 . 


2 


bisect(x 

- 3, 1, 2, 

6) 


a 

m 

b 

max_error 

1 

1.5 

2 

0.5 

1.5 

1.75 

2 

0.25 

1.5 

1.625 

1.75 

0.125 

1.625 

1.6*75 

1.75 

0.0625 

1.6*75 

1.71*75 

1.75 

0.03125 

1.71*75 

1.734375 

1.75 

0.015625 

. 1.71*75 

1.7265625 

1.734375 

0.007*125 






a + b b - a 






a , t o, 

j n 


2 2 

J 

J 


I created the pictures below working with epicycloids and hypocycloids. See 
more about this in one of the next DNLS. Josef 
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2. Window Functions 


As already mentioned, the truncation of the impulse response of a filter with a discon- 
tinuous frequency response leads to a filter whose frequency response has large oscillations in 
the vicinity of the discontinuity. Truncation of the impulse response can be considered as 
term by term multiplication of the impulse response with the rectangular window function: 

0 , 0 , 1 , ..., 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 ,..., 1 , 0 , 0 ,...} 

x 


{a n , ..., a 2 , ai, ao, a_i, a_ 2 , ..., a. n } 


or pictorially 


1 


ii 


This interpretation of truncation lends itself to extension. We can remove the oscilla- 
tions in the vicinity of discontinuities by truncating the impulse response less abruptly. For 
instance by term by term multiplication by the triangular window function: 



or a raised cosine or Hamming window: 



3. The required filters 

We will design filters to perform 3 different tasks: 

1 . A filter that will use sampled values from a function of time (sampled at unit inter- 
vals) and approximate the first derivative of the sampled function. 

2 . A filter that will use sampled values from a function of time (sampled at unit inter- 
vals) and approximate the second derivative of the sampled function. 
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#8: 


PSFTRI(f, 0, n) = MEAN(f, Q) + 


u. 


TRIW(k, n)*COFOUR(f , Q, k).EXP(- L • k *Q)1 


+ 


n 

I 

k=l 


TRIW(k , n)«COFOUR(f , Q. -k).EXP(L-k-Q) 

These functions give the frequency response of the filter found by truncation of the in- 
verse Fourier transform by the rectangular, Hamming and triangular window functions of 
length n. In other words PS FOUR represents the actual frequency response obtained by the 
FIR (= Finite Impulse Response) filter with multipliers {•••. 0, 0, a. n , •••, a n , 0, 0, •••}, and 
PSFHMG and PSFTRI represent the actual frequency responses when obtained by weighting 
the coefficients with the Hamming and triangular window functions. 

to: PIULTRECT (f , Q, n) := VECTOR(IF(k = 0, PIEAN(f , Q) , COF0UR(f, Q, k» , k, -n, n) 

#10: nULTUHCtf, Q, n) VECTOR(IF(k = 0, PIEAN(f , £1), COF0UR(f, Q, k) -HAPI(k, n», k, -n, n) 

*11: llULTTRKf, Q, n) VECTOR(IF(k = 0, FIEAH (f , £1), COF0UR(f, Q, k)-TRIW(k, n)), k, -n, n) 

These functions give the multipliers {a. n , — , a n } of the filter, using the rectangular, 
Hamming and triangular window functions respectively. 


5. The Solution Using DERIVE 

To use the frequency response method outlined above to calculate the coefficients of the 
filters which are to perform numerical first and second order differentiation, we first need to 
find the frequency response of the analytic operations of first and second order differentiation. 
The standard result from the theory of analogue control systems is that frequency response of 
a linear system may be determined by substituting s = ja> into the transfer function for the 
system. The transfer functions for first and second order differentiation are s and s 2 respec- 
tively, so the required frequency responses are ]co and -£y 2 . It should be pointed out that the 
units of co are radians per second, in contrast to the units of radians per sample given for the 
digital filters. To simplify matters, we will choose a sampling rate of one sample per second, 
in which case the two measures of frequency will be equal. 

Thus, to design a digital filter to approximate the first order derivative, Author the reqi- 
red frequency response: 

f := i-Q 

To calculate the frequency response of the filter having 6 multipliers either side of ao, 
Author and Simplify the expressions: 

PSF0UR(f, Q, 6), 

PSFTRI (f, Q, 6) and 

PSFHMG(f, Q, 6). 
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The third (low pass) filter has only its amplitude response specified - the phase re- 
sponse is of less concern in this application. If we assume a phase response which is identi- 
cally zero, the frequency response and amplitude response coincide, and the filter coefficients 
can be constructed by the same method as above. We define the desired frequency response 
using the DERIVE function STEP(t) : 

f STEP(Q+7i/4)-STEP(Q-7t/4) (5) 

We then decide on the number of coefficients in the filter and author the appropriate expres- 
sions to calculate the frequency response of the filters. For example, for ten multipiers either 
side of the central multiplier, we Author: 

PSFOUR(f , 0, 10), (6) 

PSFHMG(f , 0, 10) and (7) 

PSFTRI(f, 0, 10). (8) 

Simplify these expressions, and plot them. The required frequency response is discontinuous, 
so we would expect the Hamming window to give us the best coefficients. This should be 
verified by the plots produced. The coefficients in the Hamming filter can be calculated by 
simplifying MULTHMG(f ,Q,10). 

The values of the coefficients are: 



The values of a_io to a_i are the same as aio to ai. 


6. A DERIVE Session 

The following is a complete DERIVE session in which the coefficients for both a first 
order and second order differentiator are determined. 

We don't repeat the first 1 1 expressions which are given on pages 21 and 22 and pro- 
ceed with defining the required frequency response. 
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#28: 

#29: 


#30 

#31 

#32 

#33 


#34: 

#35: 

#36: 

#37: 

#38: 

#39: 

#40: 

#41: 


PSFHMG(f, Q, 6) 

( 2 • SIN(6*Q) f 23*73 27 31.SIN(4.Q) 9.SIN(3.Q) 

- L • I + SIN(5 *Q) + + 

l 75 l 250 125 J 200 25 

77 • SIN(2 *Q) f 23.^3 27 

+ -SIN(Q) 

100 l 50 25 J J 

Again purely imaginary, so plot imaginary part 

IM(PSFHMG(f, Q, 6)) 

and plot imaginary part of f on same axes 
IM(f) 



For a second order derivative, author appropriate frequency response 
2 

g := - Q 

nULTRECTtg, 0, 6) 

Plot and compare required and actual frequency response 

2 

12121 n 1212 1 

- j j j * j 2, — , 2, — , , — , , — 

. 18 25 8 9 2 3 2 9 8 25 IS 

PSFOURfg , Q, 6) 

2 

COS (6 ■ Q) 4 - COS (5 -Q) CQS(4-Q) 4-CQS(3-Q) n 

+ + C0S(2 -Q) + 4 ■ COS (Q) 

9 25 4 9 3 

2 

- G 







I n n 

f, Q, , , 

4 4 

PSFOURCf, Q, 10) 

PSFTRKf, Q, 10) 

PSFHnGCf, Q, 10) 

PSFOUR (black), PSFTRI (red), PSFHMG (blue) 


o.oi 
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Dependent Repeated Experiments 

Ales Kozubik Bratislava, Slovakia 


In my paper (DNL#10) I dealt with using DERIVE in teaching Probability Theory to solve prob- 
lems from independent repeated experiments. The following DERIVE code was used there: 


n ! 

pnotiCk, n) := 

#1: DIMENSION (k) 

n k ! 

n=l -i 


k n - k 

#2: ir(k, n, p) := COflBCn, k)-p ■ (1 - p) 

mlc(n s p) := 

If (n + 1) ■ p - 1 = FLOORCfn + l)-p - 1) 
tl: [fn + 1) -p - 1, (n + 1) ■ p] 

FLOORCfn + l)-p) 

lnrfa, p) := 

If LNfl - oe)/LN(l - p) = FLOORfLNfl - a)/LN(l - p)) 
#4: LNfl - oc)/LN(l - p) 

FL00RCLNC1 - a) /LNfl - p)) + 1 

l-l 

*5: af(i s p) := p-(l - p) 


#6: nrsfk, n, p) PNOflfk, n)-IF 


Dlflfk) l 

Din(p) = Din(k), n P ,0 
i=l i 


In this paper I would like to present the continuation of that programme to complete it for solving 
the problems about dependent repeated experiments and the repeated experiments when the probabil- 
ity of possible results is changed at each experiment independently on the result of the previous ex- 
periments. 

The case of dependent repeated experiments with two possible results can be modeled by the fol- 
lowing scheme: n things are in a box and k of them have a given quality and the others have not. If we 
draw m things from that box what is the probability that j of the drawn things will have the given qual- 
ity? The answer is given by the hypergeometric formula 



r n-k^ 

U'J 

\ m -j) 


(n\ 


which is realized in the function DRE ( j , m, n ) : 




COflBfk, i). COflBCn - k, n - j) 

#7: DREfj, k, n, n) := 

COflBCn, ih) 

#8: H¥PERGEOflETRIC_DENSITYCj , k, ih, n) 

HYPERGEOMETRIC DENSITY(j, Km, n) is implemented 

3600 

*9: DREO, 10, 10, 20) = 

46189 


3600 


#10: HYPERGE0flETRIC_DENSITYC3 > 10, 10, 20) = 


46189 
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#14: 


fOlHENSIONCp) 

IRCfk, p) := POLY_COEFPl II (ELENENT(p, 


i)-x + (I - ELEWENTCp, i))) s x, k 


Two worked problems: 

Problem 1. A box contains fifteen balls. Five of them are white four are red, three are blue, two 
are green and one is black. We draw at random five of them. What is the probability that 

a) two of the five balls will be white? 

b) no pair of the five balls have the same colour? 

c) the first white ball will be drawn in the fifth attempt? 


#15: 

Solution of problem 1, question a) 

The classical solutions: 

vd 

i— i 

** 

DRE(2 S 5, 5, 15) 


(5 

\ 

f 1(0 



#17: 

400 



/ 

l- 



400 


1001 

a ) 


'15) 



1001 

#1S: 

questions b) and c) 




■J 




#19: 

DRESQ1, 1, 1, 1, 1], [5, 4, 3, 2, 1]) 


5 • A 

F 3 • 

2- 

1 

40 


40 

b) - 


'15^ 



1001 

#20: 










1001 








#21: 

AFT (5 , 5, 15) 

10 

( 

) 

8 

7 5 



c - 








10 

15 

14 

13 

12 11 

#22: 










143 


Problem 2. There are four guns in the shooting gallery. Probabilities of hitting the target are 
0.1, 0.2, 0.3 and 0.4 for the single guns. One bullet is shot from each of them. What is the prob- 
ability that 

a) the target will be hit twice? 

b) the target will be hit at least twice? 

The classical solutions: 


#24: IRC (2 , [0.1, 0.2, 0.3, 0.4]) 

134 

#25: 

625 

#26: IRC (2, [0.1, 0.2, 0.3, 0.4]) + IRCC3, [0.1, 0.2, 

0.3, 0.4]) + HRC(4 , [0.1, 0.2, 0.3, 0.4]) 

643 

#27: 

2500 


0. 1 • 0.2 • 0.7 • 0.6 + 0. 1 • 0.3 • 0.8 • 0.6 + 
+0. 1 • 0.4 • 0.8 • 0.7 + 0.2 • 0.3 • 0.9 • 0.6 + 

a) +0.2 • 0.4 • 0.9 • 0.7 + 0.3 • 0.4 • 0.9 • 0.8 = 

_ 134 
625 

134 

+ 0.1- 0.2 • 0.3 • 0.6 + 

625 

b) +0. 1 • 0.2 • 0.4 • 0.7 + 0.2 • 0.3 • 0.4 • 0.9 + 

643 

+0.1 -0.2 -0.3 -0.4 = 

2500 


It should be no problem to transfer the functions to the TI-92/V 200, Josef. 
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Word Processing and 
DERIVE Revised 

J.M.M.C. Lopes & J. C. Tamames, Porto, Portugal 

In DNL#2 F. Schumm and J. Boehm discuss “Word Processing and DERIVE”. Both authors 
suggest the use of CAPTURE.COM to capture graphics and insert then in the text of a report. 

Later in DNL#6 D.A. Sjoestrand describes a method based on the capabilities of the 
Windows 3.1 Enhanced Mode. 

To insert a DERIVE full screen in the text there is an easiest way: we only have to run DERIVE 
under Windows until we reach the screen we want and then press Print Screen (PrtSc). One image of 
the screen is passed to the clipboard, and, once in WORD, we can insert this image in the text (WORD 
command Edit/Paste). We can still complete the figure (legends, etc) with Microsoft Draw, a standard 
facility in WORD. 

The capture of the screen through these methods is extremely rapid however it has the incon- 
venience of allowing only images with the original colours (at the Graphics Mode), which aren't al- 
ways the most appropriate. Of course, we can always choose the DERIVE colours as to present char- 
acters and graphics in black with a white background. But I prefer to apply PSP (Paint Shop Pro, 

JASC Inc, 1992) which has the advantage of allowing capturing an area or the full screen and adjust- 
ing the colours in the captured image. The steps to perform are: 

1 - Start Windows; open WORD and minimize then; 

2 - Open PSP and minimize then; 

3 - Start and run DERIVE under Windows until you reach the desired screen and then pass to the 

386 Enhanced Mode (ALT+ENTER): DERIVE passes running (at a low speed . . .) in a win- 
dow occupying about 2/3 of the screen. Now you can access, in the inferior part of the screen, 
the icons of WORD and PSP, making it easy to pass from one program to another 
(see Fig. 1); 

4 - Start PSP (click on the PSP icon) and then: 

Capture (Area or Full Screen) 

Make it an Icon 

(the DERIVE window reappears with the PSP cursor (a double cross); next you have to define 
the captured area with the mouse) 

Colours (Negative Image/Grey Scale/Gamma Correction 0.1) 

Edit Copy 

(one image of the desired area is copied now to the clipboard); 

5 - Close ore minimize PSP and DERIVE (ALT+ENTER to leave 386 Enhanced Mode, Quit if 

you want to leave DERIVE, else minimize DERIVE) and return to WORD (a click on the 
WORD icon). Afterwards you have to insert the contents of the clipboard in the desired posi- 
tion (Edit/Paste) 








DERIUE 
A n^ l hefTiri t h;<i I fts: 


Copyright CC) 19HH ctud 1S90 by 3i*fl Uart-iioiisH 
Honolulu* HauaiE 113ft 


Lopes & Tamames: Word Processing and DERIVE 


Mbfos eft Word 


Fig.l - The DERIVE window in 386 Enhanced Mode and the minimized icons of 
PSP and WORD. This figure is passed to the clipboard with the command 
Print Screen (with the original colours). 



1: SIN (x) 

2: SIN (3 x) 

3: 



Fig. 2 - An area screen captured with PSP, as described in the text (with the colours 
changed to black and white). 
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Ebene Algebraische und Transzendente Kurven (1) 
Die Kissoide (auch Zissoide) - The Cissoid of Diokles 

Thomas Weth, Wurzburg, Germany 


Die Theorie der ebenen Kurven ist seit alters her ein zentrales Thema der Geometrie. Durch die 
Methoden der Algebra der analytischen Geometrie und der Differentialgeometrie gelangte die Kurven- 
theorie im 18. und 19. Jahrhundert zu ihrer Bliitezeit. Durch die heute zur Verfiigung stehenden 
Moglicheiten auf geometrischem (GEOLOG, Cabri-Geometre, Felix) und algebraisch-analytischem 
Gebiet (DERIVE) gewinnen Kurven auch wieder an Bedeutung fur den Unterricht. Sie gestatten wie 
kaum ein anderes Themengebiet eine enge Verbindung von geometrischen und analytischen Methoden 
wie sie von Felix Klein zu Beginn des Jahrhunderts als „Basis“ fur den gesamten Mathematikunter- 
richt gefordert wurde. 

Da sich vor allem DERIVE dazu eignet, die Parameterdarstellung von Kurven graphisch darzu- 
stellen, die Kurven analytisch zu untersuchen, wenn moglich die algebraischen Gleichungen zu erhal- 
ten u.s.w. werden in den folgenden Nummern des DERIVE Newsletters jeweils einige der bekanntes- 
ten algebraischen und transzendenten ebenen Kurven vorgestellt. Als Literaturquelle werde ich im 
Wesentlichen die hervorragenden Bucher von Loria (Algebraische und transcendente ebene Kurven, 
Leipzig, Teubner, 1902) und von H. Schmidt (Ausgewahlte hohere Kurven, Wiesbaden, Kesselring, 
1949) zu Grunde legen (diese Werke werden im Folgenden in den Literaturverzeichnissen nicht mehr 
aufgefuhrt.) 


Die Kissoide (auch: Zissoide) - The Cissoid 


Zur Losung des Delischen Problems (Wiirfel- 
verdoppelung) (vgl. Breidenbach, 1953) wurde 
die Kissoide des Diokles (ca. 200 v. Chr.) er- 
funden. Der Name beschreibt die Ahnlichkeit 
der Kurve mit der Spitze eines Efeublattes 
(, kiggo £ Efeu). Im 17. Jhdt beschaftigten sich 
u.a. Newton, Roberval, Fermat, Huygens und 
Sluse (Domherr zu Luttich) mit dieser Kurve. 
Huygens schlug vor, die Kurve Slusianische 
Kissoide zu nennen, in der Literatur hielt sich 
aber die Verbindung des Namens zu Diokles. 


The Cissoid of Diokles (200 bC) was in- 
vented to solve the ancient problem of 
doubling a die. Its name describes the 
similarity of the curve with an ivy leaf 
( kiggo £ ivy). In the 17 th century among 
others Newton, Roberval, Fermat, Huy- 
gens and Sluse (canon in Liege) dealt with 
this curve. Huygens proposed to name the 
curve Slusian Cissoid but in literature the 
connection with Diokles was kept. 


Die uspriingliche Konstruktion der Kurvenpunkte: 

Gegeben ist ein Kreis mit dem Durchmesser |SC| = a. Senkrecht 
zu SC werden zwei zum Mittelpunkt symmetrisch liegende Seh- 
nen BBi und DDi konstruiert. Die Schnittpunkte von SB und SB! 
mit DDi ergeben zwei Kissoidenpunkte P und Pi. 

Given is a circle with diameter |SC| = a. Two chords BBi 
and DDi perpendicular to the diameter and symmetric to 
the center of the circle are drawn. The intersection points of 
SB and SBi with DDi give two points of the cissoid P and 
Pi. (See the GeoGebra realisation -» ) 
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Die Koordinaten der Kissoidenpunkte (und damit die Parameterdarstellung der Kurve) leitet man fol- 
gendermaBen her: Mit SC als x-Achse und der Kreistangente in S als y-Achse gilt fur die Hohe im 

rechtwinkligen Dreieck SBC: BF = ^x(a - x) . Mit dem Verhaltnis y : x = BF : (a - x) erhalt man: 


c(a -x) 


und daraus die algebraische Gleichung fiir die Kissoide: y 2 (a - x) - x 3 = 0 . 


Mit x = r cos cp und y = r sin cp erhalt man die Polardarstellung: r = a sin cp tancp 

Eine Parametrisierung in kartesischen Koordinaten, welche die Kurve iiber den Kreis hinaus erweitert, 
und die neben der Polardarstellung fiir DERIVE die geeignetste ist, lautet: 

at 2 at 3 

x = y = T- 

1 + t 2 * 1 + t 2 

How to derive implicit, polar and parameter form of the curve: Take SC as x-axis and the tangent of 
the circle in point S as y-axis. Then altitude BF in triangle SBC is given by BF = Jx(a-x) . Using the 


ratio x = BF : (a - x) results in y - 


:(a-x) 


which leads to the implicit form y 2 [a — x) — X 3 = 0 . 


Substituting x = rcos cp and y = r sin cp gives the polar form: r = a sin cp tancp. 


Parameterisation in Cartesian coordinates gives a parameter form: 


It might be a challenge for the students to find the given parameter form - or another one - and check 


if is “compatible” with the implicit form, (eg x = t 2 ,y 




Breidenbach, W. Das Delische Problem, Stuttgart (Teubner) 1953 

Oettinger, E. Die Zissoide oder „Efeuartige“, Mathe-Plus (Feb. 1985) Nr. 3, p. 4-7 

Weth, Th. Ein abbildungsgeometrischer Zugang zu algebraishen Kurven 3. und hoherer Ordnung, Didaiktik der 
Mathematik 19, 1991, H.2. S. 145-164 
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Gegeben: Wiirfel mit Kantenlange b und Volu- 
men b I * 3 . 

Gesucht: Konstruktion einer Kantenlange X, so 

dass V=X 3 = 2 b 3 . 

Das Problem lasst sich zuriickfuhren auf die Auf- 
losung der doppelten Proportion: 

b:X=X: Y=Y:2b , 

welche beim Auflosen A 3 = 2b 3 ergibt. Der Kis- 
soidenkonstruktion (s. oben) entnimmt man: 

Da APSDi und ASBC rechtwinkelig sind, gilt 


Cube with edge length b and 
volume b 3 . 

Finde edge length X such that 
V = X 3 = 2b 3 . 

The problem can be reduced to solving the 
double proportion: 

b : X = X : Y= Y : 2b, 

which leads to X 3 = 2b 3 . We can derive from 
the construction of the cissoide (see above): 

As APSDi and ASBC are both rigth triangles 
(with BF = DiE = z): 


(mit BF = DiE = z): 


y : x = x : z in APSDi 
x : z = z : (a- x) in ASBC 


Die Wiirfelverdoppelung mit der Kissoide 

Given: 

Searched: 


y : x = x : z = z : (a - x) 


Vergleicht man mit obiger doppelter Proportion, 
so sind PE = y = b und CE = SF = a - x = 2b als 
gegeben anzusehen. Zu konstruieren ist also die 
erste der beiden mittleren Proportionalen 
SE = x = X. 

Konstruktion: 

Zeichne eine beliebige Kossoide ( a beliebig), 
Zeichne P’ mit SP’ = a/2, 

Verbinde P’ mit C und schneide mit der Kis- 
soide - ergibt Punkt P. 

Vom Schnittpunkt P falle das Lot PE auf SC. 
Nach dem Strahlensatz gilt dann: 

PE : EC = P’S: SC =1 : 2. 

PE sollte aber gleich b sein (und damit EC = 2b). 
Deshalb muss noch eine Streckung durchgefuhrt 
werden, die PE = y auf eine Strecke der Lange b 
abbildet und gleichzeitig x im selben MaB auf die 
Lange X vergroBert (vgl. Zeichnung). Fur dieses 
(mit Hilfe der Kissoide) konstruierte X gilt nun 

X 3 = 2b 3 . 


Comparing with the double proportion from 
above we see that PE = y = b and CE = SF 

= a - x = 2b can be taken as given. We have 
ti find x with SE = x = X. 

How to do: 

Draw any cissoid (take any a), 

Draw P’ with SP’ = a/2, 

Connect P’ with C and intersect with the 
curve giving point P. 

Draw the perpendicular line PE. 

According to the theorem of proportional 
segments we finally have: 

PE : EC = P’S : SC = 1 : 2. 

PE should equal b (and EC = 2b). We need 
to find a stretching, mapping PE = y into b 
and simultaneously increasing x to X. For 
this X - found by means of the cissoid) now 
is valid: 

X 3 = 2b 3 . 


I tried to translate the important parts into English and I hope that the nongerman readers are able to 
follow Mr. Weth's interesting contribution. Josef 1993 

In 1992 Thomas Weth added a sketch of the “Delian Problem”. Now we can use DERIVE 6 and its 

slider bars to demonstrate the procedure more general together with confirming it algebraically. Josef 
2006 


p36 


Thomas Weth: Lexicon of Curves - The Cissoid 


D-N-L#ll 


* 1 : 


kissoide := 


2 3 

a ■ t a ■ t 

2 * 2 

1 + t 1 + t . 


*2: VECTORCkissoide, t s -3, 3, 0.01) 


#3: 



2 


a 


4 


*4: 

#5: 

*6: 


#7: 


S := [0 S 0], C := [a, 0] , P_ 
[P_, C] 

SOLVE (k-is son de = C + t_ ■ (P_ 

2/3 2/3 

2 2 



- C) 5 [t , t_] , Real) 

1/3 

2-2 4 


t = A t_ = 

2 


+ 

5 


Wie would like to find the solution of 
X 3 = 2 b 3 with generalized a and gener- 
alized b. 

Let's take b = 4 as an example. 


Question for the students: Why is t = 2 A (-1/3)? 


#8: P := SUBST 


2 3 

a ■ t a ■ t 


2 2 
1 + t 1 + t J 


, t, 2 


- 1/3 


#9 


: P [V T Q ’ [V "]] 


#10: [E, Q] 

#11: [S, P] 

#12: aux := Q + t ■ (S - P) 


(SOLUTIONS (a ux = 0 S t)) 

#13: 2 

1 

2977074 -b 

#14: 

913235 -a 

f 2977074 -b 

#15: E_ := SUBST aux, t, 

l 913235 -a 



[ 575302 .(913235 -a - 2977074= b) 

■ 0 

1359384087195 

#17: VECTORQE + t-(E_ - E) , 0], t, 0, 1, 0.01) 

#18: I E_ - E | 

-i4 I 13 | 

#19: 5. 502798679-10 ■ | 7 -a - 2 . 28960048- 10 -b| 

-14 I 13 I 

#20: 5.502798679-10 - | 7 -6 .4 - 2 . 28960048- 10 -4| 

#21: 5.039684198 

-!4 I 13 [ 

#22: 5.502798679-10 ■ | 7 ■ 1 - 2 . 28960048- 10 -4| 

#23: 5.039684198 

3 3 

#24: SOLVE (x = 2-4 , x, Real) 



#25: 


x = 5.039684199 
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The DERIVE Screen on the TV-Screen 

H. Scheuermann & W. Eichenauer, Hofheim, Germany 

When I was in Wolff enbiittel last year, Mr Scheuermann gave a DERIVE demonstration on a TV-screen 
using a selfmade interface. I was very much impressed and asked Mr Scheuermann to publish details 
about this interesting possibility to demonstrate PC-screens instead of using panels. Here is Mr Scheuer- 
mann 's report. I didn 7 try to translate because I am no technician, so I don 7 know the special technical 
expressions, but I am sure that you will find a way to translate the description if you are interested in 
producing this interface. Mr Scheuermann wrote about the cost: 60 German Mark. J.B. 

Zur Demonstration von Funktionsgraphen bzw. Textzeilen aus DERIVE im Rahmen von Vorfuhran- 
gen und nicht zuletzt Unterricht ist das Bild eines herkommlichen PC-Monitors mit einer Bildschirm- 
diagonalen von 14“ bzw. 16“ in der Regel zu klein. GroBbildmonitore mit Bildschirmdiagonalen von 
24“ (und mehr) geben die Bildschirmausgabe in hochster Auflosung wieder; jedoch eignen sich diese 
auf Grand ihres enormen Gewichts nur fur einen festen Standort. Dariiber hinaus ist diese Losung - 
auch heute noch - recht kostspielig. Monochrome LCD-Displays erweisen sich auf Grand der schlech- 
ten Unterscheidbarkeit von den einzelnen Graphen zB bei der Darstellung von Kurvenscharen (zB bei 
Parametervariationen) aber auch bei der grafischen Ausgabe von Funktionen als problematisch. Farbi- 
ge LC-Displays sind als Losung in der Anschafung ebenfalls kostenintensiv. Im Allgemeinen erfor- 
dern diese Displays auch noch einen speziellen Arbeitsprojektor mit geringer Warmeabgabe und eine 
gewolbte, stark reflektierende Projektionswand (haufig muss auch dann der Raum an hellen Tagen 
noch - leicht - verdunkelt werden.) Damit schlieBt auch diese Alternative einen beweglichen Einsatz 
weitgehend aus. Deshalb haben Kollegen an unserer Schule eine Schnittstelle entwickelt, die das Bild 
des Computers an den Eingang eines Farbfemsehers anpasst und iibertragt. Das Gerat ist auf Grand 
seiner kleinen Abmessungen und seiner schnellen Einsetzbarkeit (wird nur von „auBen“ an PC und TV 
angeschlossen) speziell fur den flexiblen Einsatz in verschiedenen Klassenraumen konzipiert worden. 
Da sich das Gerat auch im Eigenbau zusammensetzen lasst, ist es relativ preiswert; es verfiigt aber 
nicht fiber alle Moglichkeiten der vorher beschriebenen Alternativen. Fur viele PC-Anwendun-gen vor 
einer groBeren Anzahl von Zuschauern ist die Darstellung des Textmodus oder der Grafik in der gro- 
ben Auflosung von 600 x 200 Punkten via Adapter und handelsublichem Fernsehgerat eine zufrieden 
stellende Losung. 

Fur den interessierten Leser folgt eine kurze Beschreibung der technischen Voraussetzungen fur den Betrieb des 
Gerats und der Schaltplan der Schnittstelle *): 

PC: Die im PC eingesetzte Grafikkarte muss einen RGB- Anschluss (9-polig Sub-Min D) aufweisen. 

Bei CGA- und EGA-Karten ist dies immer der Fall. Auch an alteren VGA-Karten befmdet sich 
vielfach dieser Anschluss. Bei den neueren SVGA-Karten ist leider nur noch ein 15-poliger An- 
schluss vorhanden, sodass hier keine Anschlussmoglichkeit mehr besteht. 

TV-Gerat: Bei den meisten Geraten ist mittlerweile eine Euro-Scart-Buchse (21-polig) vorhanden. 

Software: Es konnen nur solche Programme verwendet werden, die eine Auswahl des Grafikmodus durch 

den Anwender zulassen. 

Installation: Ein TV-Gerat stellt ein Bild mit einer Zeilenfrequenz von 15625 Hz dar. VGA-Grafikkarten 

haben je nach Auflosung Zeilenfrequenzen bis iiber 40000 Hz. Um nun das Bild des PC auf ei- 
nem TV-Gerat darzustellen, muss es softwaremaBig gelingen, die Grafikkarte auf eine niedrige 
Zeilenfrequenz zu schalten. Da dies von der Grafikkarte und der verwendeten Software abhan- 
gig ist, konnen keine allgemein giiltigen Verfahren angegeben werden. Nachfolgend zwei Instal- 
lationsbeispiele: 

DERIVE 2.5: Mit Option Display EGA Color erhalt man eine Auflosung von 640x200 Punkten bei 16 Farben. 
Damit ist allerdings nur DERIVE auf dem TV-Gerat darstellbar. 

Treiberdisk.: Bei den Grafikkarten liegt meist eine Diskette des Herstellers mit verschiedenen Hilfsprogram- 
men bei. Mit einem Programm (zB CGA.COM) kann u.U. der Bildschirmmodus fur alle Pro- 
gramme umgestellt werden. 
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*) An dieser Stelle mochten wir dem Entwickler der Schnittstelle, Herrn Wolfram Eichenauer, fiir die Erlaubnis 
zur Veroffentlichung des Schaltplans danken. 





















